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Abstract. We analyze the probabilistic variance of a solution of 
Liouville's equation for curvature, given suitable bounds on the 
Gaussian curvature. The related systolic geometry was recently 
studied by Horowitz, Katz, and Katz in |12J, where we obtained 
a strengthening of Loewner's torus inequality containing a "defect 
term" , similar to Bonnesen's strengthening of the isoperimetric in- 
equality. Here the analogous isosystolic defect term depends on the 
metric and "measures" its deviation from being flat. Namely, the 
defect is the variance of the function / which appears as the confor- 
mal factor expressing the metric on the torus as f^{x, y){dx'^ +dy'^), 
in terms of the flat unit-area metric in its conformal class. A key 
tool turns out to be the computational formula for probabilistic 
variance, which is a kind of a sharpened version of the Cauchy- 
Schwartz inequality. 



Contents 

1. Liouville's equation for Gaussian curvature 

2. A geometric application 

3. Rotationally invariant case 

4. Averaging h = log / respects the differential inequality 

5. The effect of averaging on variance 

6. Some curvature estimates and Liouville's equation 
References 



1. Liouville's equation for Gaussian curvature 

Given a function / satisfying Liouville's equation for curvature in a 
domain, we are interested in studying lower bounds for the probabilistic 
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variance of /, or more precisely, for the average square deviation of / 
from its mean, in the domain. 

A geometric apphcation we have in mind is obtaining lower bounds 
for the variance of the conformal factor / in a fundamental domain of 
a doubly-periodic metric p{x,y){dx'^ + dy"^) in the plane, and hence 
for the isosystolic defect in Loewner's inequality for the corresponding 
torus (see Section [2]). 

Liouville's equation for curvature is usually stated in terms of the 
Laplace-Beltrami operator Alb on a surface with a metric p{x,y)ds'^, 
where ds^ = dx"^ + dy"^. In isothermal coordinates {x,y), the operator 

is given by A^^ = jj + ^ j . We obtain the following form of 

Liouville's equation, see [9], [131 P- 26] for details. 

Theorem 1.1 (Liouville's equation). The Gaussian curvature K = 
K{x,y) of the metric p{x,y){dx'^ + dy'^) is minus the Laplace-Beltrami 
operator of the log of the conformal factor f : 

K = -A^^\ogf. (LI) 

In terms of the flat Laplacian Ao(/i) = trace Hess (/;,), the equation 
can be written as follows: — Aglog/ = Kf"^. Setting h = log/, we 
obtain yet another equivalent form, —AqH = Ke^^. 

This equation in the case of constant curvature K is called Liouville's 
equation in [9], p. 118], see also Rogers and Schief [HI P- 154] (where it 
appears as the first Gauss-Mainardi-Codazzi equation), and [T^ [T5 | [2T]. 

Liouville's equation K = —Alb log / can be written as follows in 
terms of partial derivatives: 

- /(/.. + fyy) + Ul + fy) = K f\ (1.2) 

in other words — /Aq/ + |V/p = Kf^. The fiat Laplacian can be 
written as follows in polar coordinates: 

AoJ = ^ + + ^A^i/, [1.6) 
where A51/ is the Laplace-Beltrami operator on the circle. 



2. A GEOMETRIC APPLICATION 

In this section, we present a differential-geometric application of 
lower bounds for the variance of a solution of Liouville's equation for 
curvature. The reader mostly interested in Liouville's equation itself 
can skip to the next section. 

The systole of a compact metric space X is a metric invariant of X, 
defined to be the least length of a noncontractible loop in X. We 
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will denote it Sys = Sys(X), cf. M. Gromov [TOl E]. When X is 
a graph, the invariant is usually referred to as the girth, ever since 
W. Tutte's article [20]. C. Loewner proved his systolic inequality 
Aiea.{g) — ^Sys((y')^ > for the torus {T'^,g) in 1949, as reported 
by Pu [ig. 

The classical Bonnesen inequality from 1921 is the strengthened 
isoperimetric inequality L'^—4ttA > 7r^(i? — r)^, see [Sj p. 3], [6j. Here A 
is the area of the region bounded by a closed Jordan curve of length 
(perimeter) L in the plane, R is the circumradius of the bounded re- 
gion, and r is its inradius. The error term 7r^(i? — r)^ on the right hand 
side of Bonnesen's inequality is traditionally called the isoperimetric 
defect. 

Loewner's torus inequality can be similarly strengthened, by intro- 
ducing a "defect" term a la Bonnesen. If we use conformal represen- 
tation to express the metric g on the torus as p{dx^ + dy'^) with 
respect to a unit area flat metric on the torus R^/L (see below), then 
the defect term in question is the variance of the conformal factor / 
above. Then the inequality with the defect term looks as follows: 

Area - ^Sys^ > Var(/). 

Is there a geometrically meaningful estimate for the systolic defect term 
in Loewner's torus inequality? Thus, one could look for lower bounds 
for the variance of the conformal factor / of the metric, in terms of some 
curvature conditions, say if the curvature is bounded away from zero 
on a region D whose area is bounded below. Liouville's equation for 
Gaussian curvature K is —A log / = Kp. One is led to the following 
problem in connection with Liouville's equation. The solutions in the 
constant curvature case are of the form / = j^l'^J^'^jjp where a{z) is a 
holomorphic function on a disk (here the curvature, assumed constant, 
is normalized to the value +4) . One seeks lower bounds for the variance 
of /. Here a lower bound for the area of the region D translates into a 
lower bound for the norm of , |° IL . 

Recall that the uniformisation theorem in the genus 1 case can be 
formulated as follows: For every metric g on the 2-torus T^, there ex- 
ists a lattice L C and a positive L-periodic function f{x,y) on 
such that the torus (T'^,g) is isometric to (M.'^ /L, pds"^), where ds'^ = 
dx^ + dy"^ is the standard flat metric of M^. Similarly to the isoperi- 
metric inequality, Loewner's torus inequality relates the total area, to 
a suitable 1-dimensional invariant, namely the systole, i.e., least length 
of a noncontractible loop on the torus {T'^,g): 

Area(^7) - ^Sys((7)2 > 0. (2.1) 
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In analogy with Bonnesen's inequality, there exists a following version 
of Loewner's torus inequality with an error term: 

Area(^?)-f Sys(^7)'> Var(/), (2.2) 

see [12]. Here the error term, or isosystolic defect, is given by the 
variance Var(/) = /y2(/ ~ '^)^ of the conformal factor / of the met- 
ric g = Pqq on the torus, relative to the unit area flat metric in the 
same conformal class. Here m = / is the mean of /. More con- 
cretely, if iT^,go) = MV^, where L is a lattice of unit coarea, and D is 
a fundamental domain for the action of L on by translations, then 
the mean can be written as m = jj^f{x,y)dxdy, where dxdy is the 
standard measure of M^. 

Question 2.1. Unlike Bonnesen's inequality where the error term has 
clear geometric significance, the error term in (12. 2p is of an analytic 
nature. It would be interesting to obtain a lower bound whose geomet- 
ric significance is more transparent. Can such a bound be expressed in 
terms of suitable curvature bounds? 

The proof of inequalities with isosystolic defect relies upon the com- 
putational formula for the variance of a random variable in terms of ex- 
pected values. Keeping our differential geometric application in mind, 
we will denote the random variable /. Namely, we have the formula 

E,{f)-{E,{f)f = VaT{f), (2.3) 

where // is a probability measure. Here the variance is Var(/) = 
Ef_t ((/ — mY), where m = E^{f) is the expected value (i.e., the mean). 

Now consider a flat metric Qq of unit area on the 2-torus T^. Denote 
the associated measure by fi. Since /i is a probability measure, we can 
apply formula (12.31) to it. Consider a metric g = pg^ conformal to the 
flat one, with conformal factor f > 0, and new measure pfi. Then we 
have Efj,{p) = J^a /^/W = Area(5'). Equation (12. 3p therefore becomes 

Area((?)-(E^(/))2 = Var(/). (2.4) 

Next, one relates the expected value E^{f) to the systole of the met- 
ric g. Recall that the first successive minimum, Ai(L) is the least 
length of a nonzero vector in L. The lattice of the Eisenstein integers 
is the lattice in C spanned by the elements 1 and the sixth root of 
unity. To visualize this lattice, start with an equilateral triangle in C 
with vertices 0, 1, and | + and construct a tiling of the plane by 
repeatedly reflecting in all sides. The Eisenstein integers are by defini- 
tion the set of vertices of the resulting tiling. The maximal ratio 

o o Area 

(the area is that of a fundamental domain) for a lattice in the plane 
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is 72 = = 1.1547.... The corresponding critical lattice is homo- 
thetic to the Z-span of the cube roots of unity in C, i.e., the Eisenstein 
integers. This result yields a proof of Loewner's torus inequality for the 
metric g = Pqq using the computational formula for the variance. Let 
us analyze the expected value term E^{f) = ff^ in (12. 4p . Indeed, 
the lattice of deck transformations of the flat torus go admits a Z-basis 
similar to {t, 1} C C, where r belongs to the standard fundamental 
domain. In other words, the lattice is similar to Zr + Zl C C. Consider 
the imaginary part Q'(t) and set o"^ := Q'(r) > 0. From the geometry 
of the fundamental domain it follows that cr^ > with equality if 
and only if r is the primitive cube or sixth root of unity. Since go is 
assumed to be of unit area, the basis for its group of deck tranforma- 
tions can therefore be taken to be {a~^T,a~^}, where 3((T~^r) = a. 
With these normalisations, we see that the fiat torus is ruled by a 
pencil of horizontal closed geodesies, denoted •jy = 'jyix), each of 
length o"~^, where the "width" of the pencil equals a, i.e. the pa- 
rameter y ranges through the interval [0,cr], with = 79. By Fu- 
bini's theorem, we obtain the following lower bound for the expected 

value: E^{f) = (^J^^ f{x)dxj dy = \ength{'jy)dy > aSjs{g), see 

[T3l p. 41, 44] for details. Substituting into (12. 4p . we obtain the in- 
equality 

Area((7)-a2Sys((7)2> Var(/), (2.5) 

where / is the conformal factor of the metric g with respect to the unit 
area fiat metric go- Since cx^ > we obtain in particular Loewner's 
torus inequality with isosystolic defect (12. 2p . Hence, a metric satisfying 
the boundary case of equality in Loewner's torus inequality (12. ip is 
necessarily fiat and homothetic to the quotient of by the lattice of 
Eisenstein integers. Indeed, if a metric pds^ satisfies the boundary 
case of equality in (12. ip . then the variance of the conformal factor / 
must vanish by (12. 2p . Hence / is a constant function. 

If T is pure imaginary, i.e., the lattice L is a rectangular lattice of 
coarea 1, then the metric g = pgo satisfies the inequality 

Area((?)-Sys((?)2>Var(/). (2.6) 

Indeed, if r is pure imaginary then a > 1, and the inequality follows 
from (12. 5p . In particular, every surface of revolution satisfies (12. 6p . 
since its lattice is rectangular. 

Lower bounds for the variance of the conformal factor yield lower 
bounds for the isosystolic defect, see next section. 
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3. Rotation ALLY invariant case 

If / is rotationally invariant, then Aq/ = ^ + r§r fro'^ f!1.3p . 
Thinking of / as a single- variable function /(r, 6) = f{r), we obtain 

Ao/ = /"(r) + -/'(r). (3.1) 

r 

Liouville's equation then becomes an ordinary differential equation 

-/ (r(r) + If'ir)^ + {nr)f = Kf\r), 

cf. (11.21) . The substitution C = results in an equation in the vari- 
able (, namely 

-/T(/)+4C(/'(c)r = i^/^ 

where the relation between /(r) and f{Q needs to be explained (two 
different f's). 

Alternatively, we can proceed as follows. Instead of the substitu- 
tion above, we study the variance of the product r/(r), where / is 
the conformal factor. Thus we retain the variable r. Is there a con- 
venient form of the equation for this new function rf{r)l However, 
the natural normalisation is for J p{r)rdrd9, hence here we need 
to integrate rp{r). In other words we are integrating r/(r) with a 
weight function 1/r. We consider again the operator fl3.ip . which we 
now denote T: T = ^ + where r > 0. Applying the change 
of variable C = f"^! we can write T as T = 4^C^- Furthermore the 
substitution C = allows us to write the operator as 

4 (f 

Lower bounds for the second derivative lead easily to estimates for the 
variance, expressed by the following Lemma: If the curvature is posi- 
tive, then the function /(r) is decreasing, while log / is concave with re- 
spect to the variable t. Indeed, rewriting the equation — Tlog / = Kp 
as ~^-^2^ = ^Kf"^ immediately implies the concavity of log / with re- 
spect to t. Integrating with respect to t, we obtain — = / > 
if the curvature is positive. Hence log / is decreasing, and therefore 
so is / itself. The significance of the variable ( stems from the following 
elementary fact. 

Lemma 3.1. The variance of a rotationally invariant conformal fac- 
tor f on a disk is proportional to the variance of the corresponding 
single variable function, with respect to the variable ( = r"^. 
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Indeed, The proof is immediate from the fact that ^d^ = rdr is the 
hne measure inherited from polar coordinates. Liouville's equation for 
a rotationally invariant conformal factor / = f{C{r)) can be rewritten 
as 4^^^ log/ = —Kp. In terms of the reciprocal function = j, 

this becomes Acf^-^C,-^ log cj) = K. Note that the linear function 4>{C) = 

1 + solves this equation, as pointed out by B. Riemann, cf. f l6.ip . 

Question 3.2. Can one translate the differential inequality 

into a geometric condition involving a comparison of an arbitrary so- 
lution, with Riemann's linear solution? 

4. Averaging h = \ogf respects the differential inequality 

If the curvature satisfies a lower bound K{x, y) > a, Liouville's 
equation yields a differential inequality 

- Aq/i > ae2^("\ (4.1) 

Given a disk in where h is defined, we can average h by the 
circle of rotations about the center of the disk to obtain a rotationally 
symmetric function ha_^{r) = hav{r, 6) on the same disk, where is the 
polar angle. 



Proposition 4.1. Ifh satisfies the differential inequality ( 14. ip then its 
rotationally symmetric average h^^{r) satisfies the ordinary differential 
inequality 

- Kir) + ^hlir)) > ae'^- , (4.2) 
where a is a lower bound for Gaussian curvature. 

Proof. By linearity of the flat Laplacian, we obtain from fl4.ip — Ao(/iav) = 
— av(Ao/i) = av (_ft'e^'^) > aav (e^'^) > ae'^^''^ , by Jensen's inequality 
applied to the exponential function. Applying fl3.ip . we complete the 
proof of Proposition 14.11 □ 

The logarithmic average /la is by definition the function /la = 
gav(iog/)_ rjj^g above proposition can be restated in terms of the loga- 
rithmic average as follows. 

Corollary 4.2. // the metric f'^ds'^ admits a lower bound K > a on a 
disk, then the logarithmic average /la satisfies the differential inequal- 
ity -Alb log /la > «• 

Namely, at the level of the function h = log /, the differential inequal- 
ity —Aoh > ae^^ averages well by Jensen's inequality (see above). 
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Remark 4.3. One can find lower bounds for tlie variance of tlie con- 
formal factor /, by relating the variances Var(/) and Var(/LA), and 
developing variance estimates for a rotationally invariant function, and 
applying them to the logarithmic average /la- 



5. The effect of averaging on variance 
Lemma 5.1. Rotational averaging does not increase the variance. 

Indeed, we claim that the rotational average /lav of a function h on 
a disk D satisfies Var(/iav) < Vari^h). More specifically, let hg^^{r) = 
hUr, 0) = ^ /o'" h{r, e)de. Note that X,'" K.{r)de = 27r^ /^'^ /i(r, e)de 
lo^ /(r, 9)d6. Assume for simplicity that D has unit area. Then 

E{h^^) = j rdr (^j^ K,{r,d)dd^ = j rdr (^j^ h{r,9)d9^ = E{h). 

(5.1) 

Now consider the term EQi^). We have J^^ h1^{r)d6 = 27r/i^^(r) = 
27r /g h{r, 6)d6^ . Since the squaring function is a convex func- 
tion, we can apply Jensen's inequality to obtain (^^ J h{r,6)d6Y < 
^ J h'^{r,e)de. Thus J^"hl^{r)de < 27r^ / h'^{r,e)de. Hence 



E{hl)= I (rdr r hl{r)dd] < / {rdr r h\r)dd] = E{f' 

(5.2) 



Combining (15.11) and (15.21) . we obtain Var(/iav) = E{hl^) — E{h^-^) 
E{hl^) - E{hf < E{h'^) - E{hf, as required. 



6. Some curvature estimates and Liouville's equation 

Assuming the existence of a region of positive Gaussian curvature 
K{x,y) > a for the metric g, we would like to compare the variance 
of an arbitrary solution of Liouville's equation ( 11. ip . to that of the 
standard rotationally invariant solution /o given by 

fo = YT^2 (6-1) 

already appearing in B. Riemann's 1854 essay [T7]. 

Given an metric g = pds"^ with an arbitrary conformal factor / > 0, 
we seek estimates of the following type. 

We will denote by D{p) the disk of radius p > for the background 
flat metric of unit area. 
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Question 6.1. Suppose the torus admits a region of positive Gauss- 
ian curvature K > a, which is expressed by a conformal factor / of 
normahzed norm on D{p). We seek a lower bound for the vari- 
ance Var(/) > N{p) . . . where is an exphcit function of p. 

Applying inequalities (12.21) and f l2.5p . we can then obtain the follow- 
ing corollary: Let r be the parameter of the underlying flat metric. 
Then Area(5') — Q'(r)Sys(5')^ > N{p) .... In other words, we seek an 
estimate only dependent on the lower curvature bound and the size 
of the disk, modulo a normalisation of the L^-norm of /. Can such a 
bound be obtained by studying the reciprocal instead of /, and using 
the convexity of the reciprocal function? Is Riemann's solution fl6.ll) 
optimal as far as the variance is concerned? 

We will work with the hypothesis K{x, ?/) > a > on the met- 
ric pds^ in a disk of radius 2p > 0. We need a lower bound for the 
variance of / = /(C) with respect to the standard measure dC,, as al- 
ready discussed above in Lemma [3.11 

Lemma 6.2. Assume the Gaussian curvature K{x, y) satisfies a lower 
bound K > a > 0. Then the substitution C, = results in a con- 
cave decreasing function u{({t)) in the variable t. The second deriva- 
tive "^^2^*"*^ is bounded away from zero on the interval ( G [p, 2p] as 
follows: > ^f\ 

Indeed, since the function u = log /la is rotationally invariant, we 
can write — Tm > o/la^, or 

- ^u"{t) > «/la' (6.2) 

by (13. 2p . Monotonicity has already been checked in the Lemma above. 
Inequality (16. 2 p takes the form > jaf^. 

Can the concavity be used to obtain an estimate for the t-variance? 
Can one relate the t-variance and the (^-variance? 

References 

[1] Bangert, V; Croke, C; Ivanov, S.; Katz, M.: Filling area conjecture and oval- 
less real hyperelliptic surfaces. Geometric and Functional Analysis ( GAFA ) 
15 (2005) no. 3, 577-597. See arXiv:niatli.DG/0405583 

[2] Bangert, V; Croke, C; Ivanov, S.; Katz, M.: Boundary case of equality in 
optimal Loewner-type inequalities. Trans. Amer. Math. Soc. 359 (2007), no. 1, 
1-17. See arXiv:math.DG/0406008 

[3] Burago, D.; Ivanov, S.: Riemannian tori without conjugate points are flat. 
Geom. Fund. Anal. 4 (1994), no. 3, 259-269. 

[4] Burago, D.; Ivanov, S.: On asymptotic volume of tori. Geom. Fund. Anal. 5 
(1995), no. 5, 800-808. 



10 



MIKHAIL KATZ 



[5] Burago, D.; Ivanov S.: (Boundary rigidity and filling volume minimality of 
metrics close to a flat one. Ann. of Math. (2) 171 (2010), no. 2, 1183-1211). 

[6] Bonnesen, T.: Sur une amelioration de I'inegalite isoperimetrique du cercle et 
la demonstration d'une inegalite de Minkowski. C. R. Acad. Sci. Paris 172 
(1921), 1087-1089. 

[7] Burago, Yu.; Zalgaller, v.: Geometric Inequalities, "Nauka" (1980), Russian. 
[8] Burago, Yu.; Zalgaller, V.: Geometric Inequalities, Springer (1988). 
[9] Dubrovin, B.; Fomenko, A.; Novikov, S.: Modern geometry — methods and ap- 
plications. Part I. The geometry of surfaces, transformation groups, and fields. 
Second edition. Translated from the Russian by Robert G. Burns. Graduate 
Texts in Mathematics, 93. Springer- Verlag, New York, 1992. 
[10] Gromov, M.: Filling Riemannian manifolds. /. Diff. Geom. 18 (1983), 1-147. 
[11] Gromov, M.: Metric structures for Riemannian and non- Riemannian spaces. 

Progr. Math. 152, Birkhauser, Boston, 1999. 
[12] Horowitz, C.; Katz, Karin Usadi; Katz, M.: Loewner's torus inequality with 
isosystolic defect. Journal of Geometric Analysis 19 (2009), no. 4, 796-808. 
See arXwi0803.0690 

[13] Katz, M.: Systolic geometry and topology. With an appendix by Jake P. 
Solomon. Mathematical Surveys and Monographs, 137. American Mathemat- 
ical Society, Providence, RI, 2007. 

[14] Li, Y.; Shafrir, I.: Blow-up analysis for solutions of —Am — Ve^ in dimension 
two. Indiana Univ. Math. J. 43 (1994), no. 4, 1255-1270. 

[15] Popov, A. G.: Exact formulas for constructing solutions of the Liouville equa- 
tion A2U = e" from solutions of the Laplace equation A2V = 0. (Russian) 
Dokl. Akad. Nauk 333 (1993), no. 4, 440-441; translation in Russian Acad. 
Sci. Dokl. Math. 48 (1994), no. 3, 570-572. 

[16] Pu, P.M.: Some inequalities in certain nonorientable Riemannian manifolds. 
Pacific J. Math. 2 (1952), 55-71. 

[17] B. Riemann, Hypotheses qui servent de fondement a la geometrie. Oeuvres 
mathematiques de Riemann. J. Gabay, Sceaux, 1990. 

[18] Rogers, C.; Schief, W.: Backlund and Darboux transformations. Geometry 
and modern applications in soliton theory. Cambridge Texts in Applied Math- 
ematics. Cambridge University Press, Cambridge, 2002. 

[19] Sabitov, I.: On solutions of the equation An — /(a;,j/)e™ in some special 
cases. (Russian) Mat. Sb. 192 (2001), no. 6, 89-104; translation in Sb. Math. 
192 (2001), no. 5-6, 879-894. 

[20] Tutte, W. T.: A family of cubical graphs. Proc. Gamhridge Philos. Soc, 43 
(1947), 459-474. 

[21] Vekua, I.: Remarks on properties of solutions of the equation Am = — 2ife". 
(Russian) Sibirsk. Mat. Z. 1 (1960), 331-342. 

Department of Mathematics, Bar Ilan University, Israel 
E-mail address: katzmikfimacs . biu .ac.il 



